Vortex-Induced Phase Slip Dissipation in a Toroidal 
Bose-Einstein Condensate Flowing Through a Barrier 
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We study superfluid dissipation due to phase slips for a BEC flowing through a repulsive barrier 
inside a torus. The barrier is adiabatically raised across the annulus while the condensate flows with 
a finite quantized angular momentum. At a critical height, a vortex moves from the inner region 
and reaches the barrier to eventually circulate around the annulus. At a higher critical height, an 
anti-vortex also enters into the torus from the outer region. Both vortex and anti-vortex decrease 
the total angular momentum by leaving behind a 2n phase slip. When they collide and annihilate 
or orbit along the same loop, the condensate suffers a global 2tv phase slip, and the total angular 
momentum decreases by one quantum. In hydrodynamic regime, the instability sets in when the 
local superfluid velocity equals the sound speed inside the barrier region. 

PACS numbers: 47.32.-y,03.75.Lm,03.75.Kk 



Introduction. Flow dynamics through a constriction can 
reveal essential aspects of superfluidity. A central feature 
observed long ago with superfluid 4 He currents through 
an orifice [l[ is the occurrence of single 2n phase slips, 
which collectively decrease the fluid velocity by a quan- 
tized amount. More recently, the transition from phase 
slips to the Josephson regime has been observed by in- 
creasing the helium healing length with respect to the 
size of the orifice Q ■ 

Common belief associates phase slips with the nucle- 
ation of vortices transversally crossing the constriction 
Q . This mechanism has been invoked to explain the dis- 
sipation of the superfluid helium flow, which occurs at 
critical velocities much lower than predicted by the Lan- 
dau criterium. The microscopic mechanism of the onset 
of the instability and its dynamical evolution, however, 
are still not completely understood [J]. 

The supcrflow dynamics of a dilute Bose-Einstein con- 
densate (BEC) gas can shed new light on the physics of 
phase slips. While in quantum liquids constrictions are 
made by single or multiple orifices, in BECs they can be 
created by a laser beam generating a repulsive barrier for 
the atoms, or by an offset of the central hole in toroidal 
geometries Q. Broadly speaking, similar configurations 
allow for the observation of macroscopic phase coherence 
effects and can lead to a range of important technologies. 
While superconducting Josephson junctions are already 
employed in sensors and detectors, their superfluid coun- 
terparts can realize ultrasensitive gyroscopes to detect 
rotations Q. For instance, a toroidally shaped superfluid 
weak link provides the building block of a d.c. -SQUID, 
which is most promising sensing device based on super- 
fluid interference. 

A distinctive feature of quantum gases rests with the pos- 
sibility of experimentally interrogate the response of the 
system in a wide variety of traps and dynamical con- 
figurations. Moreover, even if the BEC is described by 
a local Gross-Pitaevskii Eq.JT]), (as in most cases where 
dipolar interactions can be neglected) and therefore lacks 



the rotonic part of the helium spectrum, its nonlinearity 
appears to be the only crucial ingredient needed to re- 
veal the microscopic mechanisms underlying the vortex- 
induced phase-slips. Superfluidity of a BEC confined in a 
torus, in absence of barriers, has been first experimentally 
observed at NIST @. The BEC was initially stirred by 
transfer of quantized orbital angular momentum from a 
Laguerre-Gaussian beam and the rotation remained sta- 
ble up to 10 seconds in the multiply connected trap. The 
metastability of a ring-shaped superflow due to centrifu- 
gal forces has been observed in Q . The superfluid critical 
velocity in a harmonically trapped BEC swept by a laser 
beam has been observed experimentally in [7( and associ- 
ated with the creation of vortex phase singularities in Q 
while solitons were observed in the one dimensional ge- 
ometry of 0] . Such problems have been object of a large 
theoretical study mainly based on numerical simulations 
of the CPE [loj]. 

In this manuscript, wc theoretically study the dynamics 
of a BEC flowing inside a toroidal trap at zero temper- 
ature and in the presence of a repulsive barrier. Similar 
qualitative results are observed when, rather than by a 
repulsive barrier, the constriction is created by an offset 
in the position of the central hole of the torus. As ini- 
tial condition, we consider a superfluid state with a finite 
orbital angular momentum in the cylindrically symmet- 
ric torus. The critical regime is reached by adiabati- 
cally raising the standing repulsive barrier. The dissipa- 
tion takes place through phase slips created by singly- 
quantized vortex lines crossing the flow. We found two 
different critical barrier heights. At the smallest critical 
height, a singly-quantized vortex moves radially along 
a straight path from the center of the torus and enters 
the annulus (Fig. QJa)-(b)), leaving behind a 2tt phase 
slip. Eventually, it keeps circulating with the background 
flow without crossing completely the torus so that it de- 
creases the total angular momentum only by a fraction of 
unity. At the highest critical height, a singly-quantized 
anti- vortex enters the torus from the outward low density 
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region of the system. The ensuing vortex dynamics de- 
pends on the velocity asymmetry between the inner and 
the outer edge of the annulus as well as on the final bar- 
rier height and ramping time. For instance, a vortex and 
an anti-vortex can just circulate on separate orbits (Fig. 
[He)) or can collide along a radial trajectory and annihi- 
late (Fig. Old)). When they orbit on the same loop or 
annihilate, the system undergoes a global 2tt phase slip, 
with the decrease of one unit of total angular momen- 
tum. In general, the BEC flow can be stabilized after 
the penetration of a few vortices. 

In hydrodynamic regime, we find that the instability to- 
wards vortex penetration occurs when the local super- 
fluid velocity equals the sound speed. This happens in- 
side the barrier region and close to the edges of the cloud. 
We have studied the above scenario in 2D and 3D numer- 
ical simulations of the dynamical GPE. In the 3D anal- 
ysis, we have employed the experimental parameters of 
the toroidal trap created at NIST [ll| . The experimental 
investigation of the system proposed here can provide the 
first direct observation of interconnection and dynamical 
evolution of vortices and phase slips in supcrfluid sys- 
tems. 

Phase-slips and vortices. We numerically solve the time 
dependent GPE 
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where g is proportional to the intcr-particle scattering 
length. In the following, we first consider an effective 
2D Cartesian geometry [12[ and eventually extend the 
analysis to the 3D configuration. The trapping potential 
Vt(r) = Vh(r) + V c (r) is made by an harmonic confine- 
ment Vh(r) = hu)±(x 2 + y 2 )/2d 2 L plus a gaussian core 
Vc(r) = Vo exp[— (x 2 + y 2 )/<r 2 ] creating a hole in the trap 
center (in what follows we will express quantities in trap 
units of time ljJ 1 and length d±). As an initial condition, 
we consider the numerical ground state obtained with 
Vb = and transfer by linear phase imprinting (in 3D 
calculations a Laguerre-Gaussian beam is implemented 
[ill ]) a total angular momentum L z = Nl, with N the 
total number of particles and I integer. The transferred 
angular momentum is low enough to have flow veloc- 
ities in the torus region much smaller than the sound 
speed. Over each loop of radius r = \J x 2 + y 2 the circu- 
lation is C = 2irl and the modulus of the fluid velocity, 
v(r) = C/2irr, is constant and directed along the tan- 
gent of the same loop. In principle, these I quanta of 
circulation can be carried by a single multiply-quantized 
macro- vortex (l3j . which hovewer breaks up into singly- 
quantized vortices still confined within the central hole 
14j . In our simulations, as soon as a finite angular 
momentum is transferred to the condensate, the vor- 
ticity field component perpendicular to the x — y plane 
v{r,t) = (V x v(r,t)) ■ z shows a "sea" of positive and 
negative vorticity spots, that is, a mesh of vortices and 
anti- vortices, Fig. [1] (b). This happens in two regions 
of very low density, close to the center and in the space 
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FIG. 1: (Color online) (a), (b) and (c) t r = 10, L z /N = 8 and 
V B — 0.34 fi. (a) t = 7.6. Density contour plot with no visible 
vortex core, (b) t — 7.6. The z component of the vorticity 
field v(r). The white dashed lines indicates the TF radii of 
the cloud. The encircled dot corresponds to a vortex about 
to enter the annulus from the inner edge, (c) t = 11.6. A 
vortex circulates along the annulus while the vorticity (inset) 
shows an anti-vortex about to enter, (d) t r = 10, L z /N = 2, 
V a = 0.61 /i. Vortex anti-vortex annihilation. 

surrounding the torus [l5j |. 

After angular momentum is transferred to the cloud, the 
barrier potential Vb(r, t) is slowly ramped up over a time 
t r to a final heigth V s . We use a repulsive well with 
widths w x centered at the maximum density and w y cen- 
tred at y = (lo| . We always choose w x > d, where 
d = R-2 — R\ is the width of the annulus. Initially, the 
density and velocity field adapt to the presence of the 
barrier, and the flow shows no sign of dissipation. In 
the barrier region, where the density is depleted, the 
flow velocity increases mainly at the edges of the an- 
nulus. By examining the vorticity, we observe that the 
two vortex seas are strongly fluctuating, with vortices 
and anti-vortices trying to escape but being pushed back 
by zones of higher density. However, when the barrier 
reaches a critical height V c \, a vortex from the inner sea 
can successfully escape and enter the annulus. As shown 
in Fig. [1] (a) and (b) [TtJ , at V c \ the flow can no longer 
sustain a stationary configuration and becomes unstable. 
In Fig. [T] (a), we observe the depletion of the density 
but not a visible vortex core. However, if we inspect the 
the vorticity field plotted in Fig. Q] (b), we clearly see an 
isolated red spot, corresponding to a positive vorticity, 
moving radially from the center of the torus towards the 
higher density region, indica ting the presence of the core 
of a singly-quantized vortex [181 ]. 

The above scenario for vortex nuclcation in a multiply 
connected geometry confirms that a persistent flow in 
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FIG. 2: (Color online) Circulation (solid lines) for loops with 
different radii and total angular momentum (dots) as a func- 
tion of time. The parameters are the same as in Fig. QJa). 
The 2ir drops in the circulation at r = 4, 4.8, 6.4 are due to 
a singly-quantized vortex moving outwards from the center. 
The drop at r = 8 and t ~ 12 is due to the passage of an 
anti-vortex entering the annulus from the outer edge. The 
oscillation in the circulation at r = 4 and t ~ 16 — 17 is due 
to a double crossing of a vortex trying to escape the inner 

region. 



such a configuration is possible because of the pinning of 
the vorticity in the low density regions near the center 
and outside of the torus. The pinning is due to the ef- 
fective energy barrier felt from a vortex core when tr yin g 
to move towards a region of much higher density [l9(. 
The effective energy barrier arises from the nonlincar- 
ity of the GPE. The obstacle raised across the annulus 
serves to unpin singly-quantized vortices by steadily de- 
creasing the density during the ramping process, up to 
suppression of the effective energy barrier. The density 
depletion occurs on a radial stripe and makes way for the 
vortex moving outwards along a straight line connecting 
the center of the vortex with the barrier. This as well 
happens for the anti-vortex moving inwards at a larger 
height of the repulsive barrier, see below. 
In the hydrodynamic regime, when £ <C d,w Xl w y and 
V s <C jU, we observe the instability towards vortex pen- 
etration when the local superfluid velocity reaches the 
sound speed [20| . This critical condition is first met in- 
side the barrier region, at the Thomas-Fermi radius of 
the cloud. The sound speed, is calculated at the maxi- 
mum of the repulsive well (at y = in our case) with the 



density integrated along the radial direction [211 ] . The 
parameters of Fig. [TJ however, have been chosen such 
that the system is outside the hydrodynamic regime, in 
order to emphasize the generality of the presented vortex 
dynamics phenomenology. 

The passage of a vortex core between two points causes 
a 27r slip in the phase difference between them 0]. In 
Fig. [2l we observe 2-7T sharp drops in the circulation C on 
a given loop of radius r at the moment the vortex core 
crosses it. Moreover, at small radii, close to the inner 
sea of vortices, the circulation shows spikes at which it 
decreases by 27t, then quickly goes back to its previous 



value. These are associated with a vortex moving out of 
the sea but being pushed back by a region of high density 
located slightly outwards, as discussed above. 
Due to phase slippage, the angular momentum is re- 
duced, and eventually the system becomes stable again 
after a finite number of spawned vortices. The circula- 
tion is lowered by a few quanta, and the fluid velocity on 
vortex-crossed loops is brought back below the critical 
value. If the ramping is stopped at V c i, only the inner 
edge of the annulus is unstable since its fluid velocity is 
larger (v(r) oc l/r) . In this case, vortices do not cross 
completely the torus and move on stable circular orbits 

However, when the barrier reaches a second critical 
height V C 2 > V c i the outer part of the annulus becomes 
also unstable. Anti-vortices then enter from outwards 
while vortices enter the inner edge, as previously dis- 
cussed. Anti-vortices move radially inwards and con- 
tribute to stabilize the outer part by phase slips. Indeed, 
an anti- vortex crossing a loop makes the circulation drop 
as a vortex crossing the opposite way. In Fig. [ljc) we 
see a vortex already circulating inside the high density 
region of the annulus while an anti-vortex begins to en- 
ter. The separation between V c i and V C 2 is proportional 
to the velocity difference Av = 1{R\ — R2)/{R\R2) be- 
tween the two edges. In general, depending on Av, the 
dynamics at barrier heights larger than V C 2 can vary. For 
instance, at lower angular momenta Aw becomes smaller, 
and a vortex and an anti- vortex enter the annulus almost 
simoultaneously. They can then collide and annihilate, as 
shown in Fig.[TJd). When a vortex and an anti- vortex an- 
nihilate or separately orbit on the same loop, the system 
undergoes a global 2ir phase slip, and the total angular 
momentum is decreased by one unit. 
We extend our 2D calculations into a 3D configuration 
[23l ] . The parameters of the toroidal trap are those em- 
ployed experimentally at NIST [ll[ . We add a repulsive 
well potential [l6[ whose shape, however, is not crucial 
in determining the qualitative features of the dissipation 
as long as w x is larger than the width of the annulus. 
Since the healing length is of the order of the harmonic 
lenght along z we found, as expected, that the nucleation 
of singly-quantized vortex lines and their dynamics re- 
semble those observed in 2D calculations. In particular, 
we have two critical values for the barrier height V c \ and 
V C 2 connected respectively with the nucleation of vortices 
or both vortices and anti- vortices. 

Conclusions. We have studied the superfluid dynamics of 
a dilute Bose-Einstein condensate confined in a toroidal 
trap in presence of a repulsive barrier. With a finite ini- 
tial angular momentum, we observed two critical values 
of the barrier height for the onset of phase slips dissi- 
pation: a lower one corresponding to vortices entering 
the annulus from the center of the torus, and a higher 
one related to both vortices and anti-vortices, the latter 
entering the outer edge of the annulus. We have per- 
formed 3D simulations with the NIST toroidal trap pa- 
rameters, where the above scenario could be experimen- 
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tally observed when a standing repulsive barrier is raised 
across the BEC supenuid flow. Since supcrcurrents have 
recently been observed in absence of the barrier, we be- 
lieve that the experimental confirmation of our results is 
at hand. Vortices can be directly observed with BECs, 
and it is therefore possible to experimentally characterize 
their role in phase-slips-induced dissipation in superfluid 
systems. 
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